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Abstract. We propose a universal approach to a range of enu- 
meration problems in graphs. The key point is in contracting suit- 
ably chosen symmetric tensors placed at the vertices of a graph 
along the edges. In particular, this leads to an algorithm that 
counts the number of G?-regular subgraphs of an arbitrary graph in- 
cluding the number of d-factors (previously we considered the case 
d = 2 with a special emphasis on the enumeration of Hamiltonian 
cycles; cf. math. CO/0403339). We briefly discuss the problem of 
the computational complexity of this algorithm. 



1. Introduction 

A tensor network is a collection of tensors placed at the vertices of 
a connected finite graph. The valency of each tensor must coincide 
with the degree of the corresponding vertex, and tensor indices are 
labelled by the half-edges incident to this vertex. Two tensors placed at 
adjacent vertices can be contracted over the pair of indices labelling the 
two half-edges of the connecting edge. This yields a tensor network with 
one tensor less than the original one (the underlying graph obtained by 
contracting the corresponding edge has also one vertex and one edge 
less). When contracting a loop, the tensor at its single vertex gets 
contracted over a pair of indices, and the total number of tensors (or 
vertices) remains unchanged. The total contraction along all edges of 
the graph gives a number (the graph reduces then to a single vertex 
with no edges). 

The idea of placing tensors at vertices of graphs and contracting them 
along the edges traces back to Sylvester [H] 1 , who used this construction 
for obtaining polynomial invariants of symmetric tensors. Much later 
Penrose |3] reversed this idea to computing graph invariants (like the 
number of 3-edge colorings of a planar 3-regular graph) in his novel 
(though not quite successful) approach to the 4-color problem. Later 
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a number of other remarkable applications was found, first of all in 
the theory of Vassiliev knot invariants (see, e. g., pQ). Recently tensor 
network contraction was used to simulate quantum computation (cf. 
[2] and references therein). 

Here we deal with contracting networks of specially chosen symmetric 
tensors. The result of contraction gives a generating function for the 
numbers of regular subgraphs of any type. In its simplest version this 
algorithm produces the number of factors of any given degree. We also 
comment on the computational complexity of this agorithm. 

Acknowledgements. We are grateful to S. Cook and M. Braver- 
man for useful remarks. 

2. Networks of symmetric tensors and graph functions 

Let G be a finite graph (possibly with loops and multiple edges). 
The set of vertices of G we denote by V(G) = {i>i, . . . , v n }, where 
n = |V(G)| is the total number of vertices, and the set of edges of G 
we denote by E(G). For each vertex Vi we denote by d\ its degree (or 
valency), i — 1, . . . , n. Then the number of edges of G is given by 

1 n 

\E(G)\ = - J £d i . 

i=l 

Now let F be a field, and let V = F r be a vector space of dimension 
r over F. Fix a symmetric bilinear form B : V <S> V — > F. The graph 
G together with the bilinear form B define a multilinear form 

(1) B G : V® dl ® ...®V® dn — >F, 

which is constructed as follows. At each vertex Vi of G we place di-th 
tensor power V® dl of the vector space V, where the factors are labeled 
by the half-edges of G incident to Vi. Each edge of G defines a contrac- 
tion of two copies of V (corresponding to its two half-edges) by means 
of the bilinear form B. We obtain the multilinear form Bq by perform- 
ing such contractions over the set E{G) of all edges of G. Rigorously 
speaking, the multilinear form Bq depends on the order of half-edges 
at each vertex v^, or, equvalently, on the order of factors in the tensor 
power V® di . However, its restriction to S dl V <E> • • • <E> S dn V, where S d V 
denotes the d-th symmetric power of V, is defined uniquely. 

Now fix a sequence A = {A\, A2, ■ ■ ■} of symmetric contravariant 
(i-valent tensors Ad £ S d V C V® d . Here we consider tensor networks 
given by the triple {G, A, B}. We treat the tensor product A^ ® • • • <g) 
Ad n as an element of V® dl <S> ■ ■ ■ <S> V® dn and consider the element 

(2) Fa,b(G) = B G (A dl ® • • • ® A di ) e F. 
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In other words, J-^b(G) is the result of contracting the tensor network 
{G, A, B} - it is obtained by placing a copy of A d at each vertex of 
G of degree d and contracting ®™ =1 A d . using B over \E(G)\ pairs of 
indices corresponding to the edges of G. Thus, to each pair A, B, 
where A is a sequence of symmetric d-tensors (d — 1,2,...) and B is 
a symmetric bilinear form, we associate an F- valued mapping JF4 B on 
the set of isomorphism classes of graphs, or an F- valued graph function 
in the terminology of jHj. 



Given a graph G, let H be a ^-regular subgraph of G (not necessarily 
connected). The type of the subgraph H in G is the partition A# = 
[|y(#i)|, . . . , \V(Hi)\] of the number \V(H)\ of vertices of H, where 
Hi, . . . , Hi are the connected componets of / = 1(H). The weight 
of partition A# is |A#| = and the length is 1(Xh) = 1(H)- For 

each partition A we define a graph function by 

Nk,\(G) = #{H C G|# is ^-regular and A^ = A}, 

i.e., Nk,\(G) is the number of ^-regular subgraphs of type A in G. 
Clearly, Nk,\(G) — if k > max{di, . . . , d n } or |A| > n. 

For k a positive integer, denote by Pk(%i, x 2 , ■ ■ ■) = x\ + x\ + . . . 
the fc-th power sum in variables Xi,X2,.... Given a partition A = 
[ki, . . . , ki\, we define a homogeneous symmetric function p\ of degree 
I A I = ki + ■ • • + k[ by the formula 



We want to show that under a special choice of A and B the graph 
function T^b defined in Section |2 counts the number of fc-regular 
subgraphs of any given type in graphs. We take F = C and consider 
the standard coordinates in V = C r . In these coordinates the bilinear 
form B is given by the identity r x r matrix I r . We define the tensors 
A^k componentwise. For d < k we put 



3. Enumeration of regular subgraphs 



p x (x u x 2 ,...) = Y[ Pki(xi,x 2 , . . .). 



i=l 
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Ah — id — . 
A dM - N 



and for d > k we put 

Xi if (ii . . . id) is a permutation of r_^r), 

k d—k 

% = l,...,r - 1, 

* if (h . . . i d ) = (r . . . r), 
otherwise, 

where Xi,...,x r _i and t are arbitrary complex numbers. The main 
result of this section is the following 

Theorem 1. For Ak = {^4i,fe, ^-2,fc, • • • } as above, the value of the graph 
function J-A k ,i r on an V graph G is given by the formula 

(3) FA k ,i r (G)= ^r^i.-.^^lG), 

|A|<n 

where the sum is taken over the set of all partitions A of weight \X\ < 
n=\V(G)\. 

Proof. We interprete the indices 1, . . . , r as colors of the half-edges of G. 
A product of n components ^ d \ k %dl ■ ■ ■ ^dZ~k" +1 '" lm ( wnere n — \V(G)\ 
and m = YTj=i dj — 2\E(G)\) contributes non-trivially to J r A k ,i r ( < ^) if 
and only if the colors agree on each edge of G or, equivalently, if and 
only if for each edge the both indices that label two of its half-edges 
are the same. Thus, in this case the non-zero contributions are in one- 
to-one correspondence with edge colorings of G in r colors with the 
following properties: 

(i) an edge incident to a vertex of degree < k has color r, and 

(ii) at each vertex Vj G V(G) of degree dj > k exactly k edges 
incident to it have some color ij G {1, . . . , r}, and the rest di — k edges 
have color r (if an edge makes a loop we count it twice). 

The closure of the union of edges with colors 1, . . . , r — 1 is a k- 
regular subgraph H in G, and every connected component Hj of H 
is colored in one of the colors ij G {1, ... ,r — 1}. The contribution 

to J 7 A k ,i r (G) fr° m this coloring is fH Aff l Y[ l j=i X ^ H ^\ where Xh = 
[\V(Hi)\, . . . , |V(i?j)|] is the partition associated with H and I = 1(H) is 
the number of connected components of H. Therefore, the contribution 
from all possible colorings of the subgraph H is equal to 

fHAffITT | y^x[ V[Hj}l ) =t n - lXHl px H (x 1 ,...,x r - 1 ), 




and summig up the contributions from all A;-regular subgraphs in G we 
get the assertion of the theorem. □ 
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Corollary 1. The graph function J-^. k j r , depending on xi,...,x r _i 
andt as parameters, determines the numbers N^x{G) uniquely for any 
graph G with n < r — 1 vertices. 

Proof. By Theorem ^ the graph function T_A k j r with values in C fac- 
tors through the ring C[xi, . . . ,x r -i] Sr ~ 1 of symmetric polynomials in 
r — 1 independent variables x±, . . . , x r _i. It is well known that the 
polynomials Pk(xi, • • • , 2V-i)> k = 1, ■ ■ ■ ,n, are algebraically indepen- 
dent in C[x±, . . . , x. r _i] 5r_1 provided n < r — 1. Therefore, in this case 
the graph function J-A k ,i r {G) determines the coefficients Nk,\{G) in Q 
uniquely. □ 

Remark 1. Since the coefficients Nk t \(G) in J3J) are non-negative inte- 
gers, we can uniquely find them out when 



where is the number of connected componenets of H and the 

maximum is taken over all fc-regular subgraphs H in G, but we will 
not dwell on this issue here. 

Below are two special cases of Theorem^ of independent interest. 

Corollary 2. Put r — 2, X\ — 1 and t = 0. Then for any graph G 
the value J r ^ k j 2 (G) is the number of fc-factors, or spanning k-regular 
subgraphs in G. 

Proof. By Theorem Q 



In the simplest case when k — 1, Corollary 1 counts the number of 
1-factors, or perfect mathcings in G. 

The next statement concerns connected ^-factors of G: 



Corollary 3. Put r = n + 1, Xj = e 2lTyr ~^^ n , (j = 1, . . . , n) and t = 0. 
Then for any graph G the number of connected k-f actors in G is equal 
t° ^^A k ,i n+ i{G) , where n is the number of vertices of G. 

Proof. In this case 



r > 1 + max 

HcG 



□ 




0, k — 1, . . . , n — 1 



By Theorem^ we have 
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where [n] denotes the 1-element partition of n = \V(G)\ (by definition, 
Nk,[n](G) is the number of connected fc-factors). □ 



4. Computational complexity 

The value J r A k ,i r (^) can be effectively computed for any graph G as 
explained in Section In general, this is a hard computational prob- 
lem. To make things less complicated let us take r = 2. Then Corollary 
121 provides a relatively simple algorithm for counting the number of k- 
factors in graphs. (Note that even the case k — 1, or enumeration of 
perfect matchings, is a well-known #P-complete problem.) 

The computational complexity of the above algorithm depends on 
two main points: 

(1) the succession of tensor contractions along the edges of G (or 
simply edge contractions); 

(2) the utilization of the special form of tensors A^k] 

What concerns the first point, the objective is to choose a sequence 
of edge contractions that keeps the maximal tensor valency (or vertex 
degree) as small as possible. This problem was solved in |2J. Take a 
sequence of edge contractions reducing G to a point. The complexity 
of this sequence is the maximum vertex degree during the contraction 
process. The contraction complexity cc(G) is the minimal complexity 
over all sequences of edge contractions. A nice result of j2] states that 
cc(G) = tw(G*) + 1, where G* is the linear graph of G, and tw(G*) 
is the treewidth of G*. Though computing the treewidth of a general 
graph is iVP-hard, a tree decomposition of G* of width 0(tw(G*)) can 
be obtained in time \V(G*)\° W exp(0(tw(G*)) g]. Moreover, given a 
tree decomposition of G* of width w, a sequence of edge contractions 
in G of complexity not greater than w + 1 can be found in polynomial 
time. Thus, the results of [2] give a rough estimate of the computa- 
tional complexity of our algorithm by \V(G)\°^ exp(0(tw(G*))). So, 
for graphs whose line graphs have bounded treewidth our algorithm 
works quite fast. Moreover, known constructions of almost optimal 
tree-decompositions (cf. gj) give almost optimal edge contraction se- 
quences. 

As for the second point mentioned above, the situation is less clear 
at the moment. Note that the estimate of |5j is rather general and 
applies to arbitrary tensor networks. In our case it can be considerably 
sharper, since the tensor networks we use are of very special types. In 
particular, Sylvester's theory of decomposing symmetric tensors into 
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sums of symmetric powers of vectors may appear quite useful in im- 
proving the efficiency of our algorithm. 2 However, this is still a work 
in progress. 
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